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ABSTRACT 


This  report  deals  with  the  optimum  structural  design  of  circular  cylin¬ 
drical  shells  reinforced  with  identical  uniformly  spaced  T-ring  stiffeners, 
and  subjected  to  external  pressure  loading.  The  optimization  problems  con¬ 
sidered  are  of  three  types:  (1)  minimum-weight  design,  (2)  design  for  maximum 
separation  of  the  lowest  two  natural  frequencies,  and  (3)  design  for  maximum 
separation  of  the  lowest  two  natural  frequencies  which  have  primarily  axial 
content.  Gross  buckling  is  precluded  by  specifying  a  minimum  natural  fre¬ 
quency,  and  additional  behavioral  constraints  preclude  yielding  or  buckling  of 
panels,  T-ring  stiffeners,  and  web  and  flange  instabilities  within  each  T-ring. 
Analysis  is  based  on  use  of  an  equivalent  orthotropic  shell  model,  and  optimiza¬ 
tion  is  accomplished  through  use  of  a  sequential  unconstrained  minimization 
technique.  Examples  indicate  that  a  small  increase  in  weight  above  optimum 
(minimum)  values  can  result  in  relatively  large  increases  in  frequency 
separation,  and  that  maximum  frequency  separation  is  obtained  when  second 
and  third  lowest  frequencies  approach  each  other. 
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Thickness  of  shell  skin  (in) 

Thickness  of  frame  web  (in) 

Thickness  of  frame  flange  (in) 

Depth  of  frame  web  (in) 

Width  of  frame  flange  (in) 

Spacing  of  frames  (in) 
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Mass  density  of  structural  material  (slug/in  ) 
Total  length  of  structure  (in) 

Radius  to  mid-surface  of  skin  (in) 

Hydrostatic  pressure  (psi) 

Kinetic  energy 
Shell  coordinate  system 
Midsurface  displacements 
Displacements  of  an  arbitrary  point 
Eigenvector  of  ith  natural  frequency 
Natural  frequency 
Wave  numbers 
Time  (sec) 

Stiffness  matrix 
Geometric  stiffness  matrix 
Mass  matrix 
Yield  stress  (psi) 

Critical  flange  buckling  stress  (psi) 

Critical  web  buckling  stress  (psi) 

Objective  function 
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Vector  of  design  variables 
Design  constraint 
Number  of  design  variables 

Critical  pressure  for  buckling  of  skin  (psi)  between  frames 
Minimum  allowable  natural  frequency  of  vibration  in  vacuo  (Hz) 
Young’s  modulus 
Poisson's  ratio 

Unsupported  length  of  shell  plating  (in) 

Maximum  allowable  normalized  structural  mass 
Composite  objective  function 

Constant  for  use  in  computing  extended  penalty  function 
Positive  scalar  quantity 

Direction  vector  for  uni-directional  search 
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INTRODUCTION 


Although  a  wealth  of  literature  exists  for  the  static,  dynamic  and  sta¬ 
bility  analyses  of  stiffened  shells  of  revolution  subjected  to  various  applied 
loads,  with  the  majority  of  these  studies  devoted  to  cylindrical  shells,  the 
work  of  Schmit  and  Morrow  [1]  serves  as  a  pioneering  effort  toward  the  intro¬ 
duction  of  structural  optimization  concepts  into  the  design  of  stiffened  cylin¬ 
drical  shells.  In  this  reference  a  cylindrical  shell,  reinforced  with  longi¬ 
tudinal  and  ring  stiffeners,  each  with  rectangular  cross  section,  was  designed 
to  carry  a  number  of  independently  applied  sets  of  static  loads  with  minimum 
structural  weight.  The  shell  was  constrained  against  overall  (system)  buckling, 
panel  and  stiffener  buckling,  and  also  against  material  yield.  The  mathematical 
model  which  formed  the  basis  for  the  stress  and  buckling  analyses  was  an  equiv¬ 
alent  homogeneous  orthotropic  shell;  i.e. ,  the  discrete  stiffeners  and  skin 
stiffness  properties  were  incorporated  in  the  orthotropic  elastic  shell  stiff¬ 
ness  properties.  This  theory,  a  3rd  order  Flugge-Lur 1  e-Byrne  type  theory, 
proved  adequate  provided  the  stiffener  spacing  and  cross  sectional  dimensions 
were  sufficiently  small  to  permit  the  smoothing  operation  inherent  in  the 
orthotropic  shell  model. 

In  a  more  recent  study,  Pappas  and  Allentuch  [2,3]  investigated  the  mini¬ 
mum-weight  design  of  ring  stiffened  cylindrical  shells,  subjected  to  a  number 
of  static  applied  load  conditions.  In  this  study,  the  ring  stiffeners  were 
T-shaped  rather  than  rectangular.  The  structural  analyses  of  general  instabil¬ 
ity,  localized  panel  instability,  and  stiffener  instabilities  were  based  on 
buckling  formulas  contained  in  Ref.  [4],  along  with  the  appropriate  stress 
limits . 

The  subject  of  this  report  is,  in  essence,  a  structural  optimization 
study  which  employs  both  a  combination  and  an  extension  of  the  structural 


1 


models  in  Refs.  [1-3].  In  the  present  study,  three  somewhat  different  structural 
optimization  problems  have  been  formulated.  In  the  first,  the  minimum-weight 
design  of  a  ring  stiffened  cylindrical  shell  subjected  to  several  applied  static 
loads  is  considered.  The  structure  is  constrained  against  all  of  the 
buckling  modes  considered  in  Refs.  [2,3]  although  for  some  of  the  modes,  the 
constraint  equations  differ.  Further,  the  structure  is  constrained  against 
in  vacuo  natural  vibrations  below  a  specified  frequency  limit.  In  the  second 
formulation  the  shell  is  designed  to  maximize  the  separation  between  the  lowest 
two  in  vacuo  natural  frequencies,  while  being  constrained  against  buckling  and 
yield  behavior  as  specified  in  the  first  design  problem,  and  while  having  a 
weight  less  than  a  prescribed  maximum.  The  third  formulation  is  similar  to  the 
second,  with  the  distinction  that  the  frequencies  being  separated  are  the  two 
lowest  which  have  primarily  axial  content. 

This  optimization  study,  in  either  the  weight  minimization  form  or  the 
frequency  separation  forms,  represents  a  considerable  advance  beyond  previous 
characterizations  of  the  optimization  problem.  The  structural  model  in  this 
study  is  based  on  the  equivalent  orthotropic  shell  model  of  Ref.  [1],  with 
dynamic  effects  added.  This  representation,  although  somewhat  imprecise  in  its 
ability  to  model  a  ring  stiffened  shell,  has  proved  invaluable  in  providing  an 
economical  and  yet  reasonably  complete  initial  structural  representation  for 
use  in  the  optimization  studies  performed.  A  more  detailed  description  of  the 
mathematical  model,  and  a  summary  of  the  results  of  this  research,  are  presented 
in  the  following  sections. 
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ANALYSIS  OF  STIFFENED  SHELLS 


The  investigation  detailed  in  this  report  is  concerned  with  the  optimum 
structural  design  of  circular  cylindrical  shells  with  uniformly  spaced  T-ring 
stiffeners  (Fig.  1)  and  subjected  to  several  different  applied  loading  conditions 
typical  of  submerged  vessels,  namely  (1)  specified  external  pressure  (or  vessel 
depth),  (2)  specified  axial  compressive  loadings,  and  (3)  applied  static  loads 
associated  with  vessel  motion.  Design  variables,  shown  in  Fig.  1,  are  skin 
thickness  (tg) ,  stiffener  web  thickness  (t^) ,  web  depth  (d^) ,  spacing  (Z^) ,  flange 
width  (df ) ,  and  flange  thickness  (t^) .  Radius  (R) ,  length  (L) ,  and  the  material 
properties  are  assumed  to  be  preassigned  parameters. 

All  of  the  research  performed  to  date  has  employed  the  simplified  ortho¬ 
tropic  shell  model  given  in  [1]  in  the  calculation  of  the  natural  frequencies 
of  vibration  and  in  the  overall  (system)  buckling  analysis.  This  idealization 

has  of.fered  the  advantage  of  mathematical  simplicity  and  the  disadvantage  of  a 

•k 

somewhat  limited  modeling  capability  ,  but  it  has  provided  an  economical  initial 
structural  representation  for  use  in  the  optimization  studies. 

In  order  to  analyze  the  dynamic  response  of  the  cylinder  in  Fig.  1,  it  is 
hypothesized  (as  in  Ref.  1)  that  the  frames  and  skin  act  as  a  unit  accurding 
to  the  Bernoulli-Euler  deformation  assumption  as  extended  through  the  Flugge- 
Lur ’ e-Byrne  theory,  and  that  the  stiffness  and  inertia  properties  of  the  frames 
are  uniformly  distributed  over  the  length  of  the  cylinder.  It  is  then  possible 
to  express  the  kinetic  energy  of  the  shell  in  the  form 


T  = 


L  2  IT 

f/fi 


o  12  12 

(u  +v  +w  )dzdcpdx 


(1) 


The  model  gives  a  very  accurate  representation  of  the  structural  behavior  pro¬ 
vided  the  characteristic  wavelengths  of  the  modes  of  vibration  (or  of  the  static 
displacements)  are  very  long  compared  to  both  the  ring  spacing  and  ring  cross- 
sectional  dimensions. 
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(a)  HULL  SEGMENT  CROSS-SECTION 


Figure  1.  Typical  Hull  Segment  Cross  Section  (from  Ref.  2). 
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where  u,  v  and  w  are  the  displacements  of  an  arbitrary  point  in  the  structure  in 

the  x,  <j>,  and  z  directions,  respectively  (Fig.  2),  and  are  given  by  [1J 
9w 


u  =  u-z 


3x 


z.  z  3w 
v  =  (1TT)V-  r  3? 


(2a) 
(2b) 

w  =  w  (2c) 

where  u,  v  and  w  are  the  displacements  of  the  shell’s  mid-surface.  From  the 
kinetic  energy  the  inertia  terms  in  the  appropriate  equations  of  motion  are 
obtained  by  use  of  Hamilton’s  principle.  These  terms  are  appended  to  the 
equations  of  static  equilibrium  for  the  stiffened  shell  in  [1] ,  which  also  con¬ 
tain  the  influence  of  destabilizing  forces. 

Assuming  that  the  external  loads  give  rise  to  circumferential  and  longi¬ 
tudinal  compressive  forces  per  unit  length  of  magnitude  PR  and  PR/2,  respec¬ 
tively,  where  P  is  hydrostatic  pressure,  then  the  combination  of  inertia,  static, 
and  destabilizing  forces  leads  to  the  following  three  coupled  partial  differen¬ 
tial  equations  of  motion. 


+  h*  -  -  r»"  -  px 

h +  -  ik*  -  Vs  -  k**  +  "*>  -  pv"  =  P, 


x<|>  R  x<J)  2  9  R 
K 


mx  +  hi + ht + ¥$* +  h +  =  P; 

K 


where 


P  =  P 
x  s 


j*  [ii  -  zw 


T](l  -  f)dz 


y  s 


p  =  p 
z  ’s 


P„  =  P „/[(!-  |)2v  -  (1  -  f)(f)w*](l  -  |)dz 


R 


R'  R 


(3a) 

(3b) 

(3c) 

(4a) 

(4b) 


f  [W  +  zu’  +  (1  -  |)  (|)v*  -  (z)2wm  -  (|)  w**] (1  -  |)dz  (4c) 


R  VR' 


R 
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•  g/  \  g /  \  d(  )  ! 

and  where  (  )  denotes  -Jp  ,  (  )’  denotes  ,  (  )*  denotes  -gj-  and  J  (  )dz 

is  the  integral  through  the  thickness  of  the  shell  and  frame.  The  forces  M  and 

N  may  be  expressed  in  terms  of  the  mid-surface  displacements  u,  v  and  w  (see 

Appendix  I)  so  that  Eqs.  (3)  can  be  expressed  in  terms  of  displacements  only. 

Under  the  assumption  that  the  boundary  conditions  are  of  the  simple  support 

type,  then  the  solution  to  these  equations  takes  the  form 
.  ,  mUx 

u  =  A-^sxn  n<p  cox  — j—  sin  Ci)T 

,  .  mTTx  .  ,cs 

v  =  A  ;os  nq>  sin  — sin  WT  (5) 

w  =  A^sin  n4>  sin  sin  U)T 

where  n  =  0,1,2,...  and  m  =  1,2,.... 

Substitution  of  Eqs.  (5)  into  Eqs.  (3)  gives  the  algebraic  eigenvalue  problem 
([K]  -  P[Kg]  -  i*)2[M]){A}  =  {0}  (6) 

where  the  stiffness,  "geometric"  stiffness,  and  mass  matrices  are  [K] ,  [K^J  and 
[M] ,  respectively,  and  are  given  in  Appendix  I. 

It  should  be  noted  that  the  sine  and  cosine  dependencies  on  the  angle  (p, 
and  the  similar  dependencies  on  the  axial  variable  x,  could  have  been  inter¬ 
changed  without  influencing  the  matrices  in  Eq.  (6)  for  n  >  0.  The  n  =  0  case 
as  given  in  Eq.  (6)  is  actually  a  combination  of  the  solution  form  in  Eqs.  (5) 
(pure  torsion)  and  the  similar  form  with  sine  and  cosine  terms  (with  argument 
zero)  interchanged  (torsionless  motion) . 

Since  the  algebraic  eigenvalue  problem  for  given  values  of  m,  n  and  P  is  of 
only  rank  three,  its  solution  for  the  natural  frequencies  (eigenvalues)  and 
associated  modes  (eigenvectors)  is  accomplished  without  difficulty. 

In  the  portion  of  the  structural  optimization  research  wherein  the  structure 
was  designed  for  the  greatest  separation  of  the  lowest  two  axial-type  vibratory 
modes,  the  modes  with  A^  =  1,  and  A£,  A^  <  1  were  termed  "axial."  In  order  to 
prevent  any  general  buckling  from  occurring,  all  the  frequencies  associated  with 
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values  n  =  0,...,6,  m  =  were  retained  and  forced  to  exceed  a  prescribed 

minimum  CO  .  . 

mm 

In  addition  to  this  gross  buckling,  it  is  necessary  to  be  able  to  detect 
the  occurrence  of  several  additional  modes  of  "local"  failure.  These  local 
failure  modes  consist  of  panel  (inter-ring)  yielding  or  buckling  and  yielding 
or  buckling  of  the  webs  and/or  flanges  of  the  stiffener  rings. 

Panel  (skin)  yielding  will  be  precluded,  according  to  the  Von  Mises 
criterion,  provided 


2 

4> 


-  CT.G  + 
<px 


a2)1/2 

X 


<  a 
“  y 


(7) 


where  G^  =  material  yield  stress  in  uniaxial  loading,  and  G^  and  G^  are  in¬ 
plane  stresses  normal  to  the  surfaces  of  the  element  in  Fig.  2.  From  Ref.  [4] 
the  maximum  bending  stresses  in  the  panels  due  to  external  pressure  are 


0  .  -fV  +  r(Hn  +  vhe)] 

S 

°x  ■  ~r  <^2  +  raE> 

s 


(8a) 

(8b) 


where  T,  Hg,  and  Hn  are  load  factors  defined  in  Appendix  II. 

A  suitable  approximate  formula  [4]  for  the  critical  external  pressure 
which  will  cause  panel  buckling  is 

Pcr  =  2.42E(ts/D)5/2/  {(1-v2)3/4[(H/D)  -  0.45(ts/D)1/2]}  (9) 

where  &  =  &x  -  t^  is  the  unsupported  length  of  a  panel,  and  D  =  2R. 

Again  following  Ref.  [4],  the  maximum  compressive  stress  in  the  rings  may 
be  taken  as 


Gc  =  -PRQ/A  (10) 

where  A  =  tgt^  +  t^d^  +  b^d^,  and  Q  is  a  load  factor  defined  in  Appendix  II. 

The  magnitude  of  G  must  be  less  than  the  yield  stress,  G  ,  and  also  less  than 
C  y 

the  critical  values  of  compressive  stresses  at  which  the  flange  or  web  will 
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buckle.  Assuming  that  the  web  and  flange  are  infinitely  long  rectangular 
plates,  that  the  web  is  simply  supported  along  all  edges,  and  that  the  flange  is 
simply  supported  along  three  sides  and  free  on  one  edge  (all  conservative  assump¬ 
tions)  ,  then  the  critical  stresses  for  buckling  of  the  flange  and  web,  respec¬ 
tively,  are  [5] 


0.506TT  Er^  ,f3  _  N12 

'f  =  - 2~[2t  /(d  -t  )] 

cr  12 (1-v  )  v 


(11a) 


o 

w 


cr 


/  2_ 

4tt  E 

12 (1-V2) 


[W 


(Hb) 


It  should  be  noted  that  Eq.  (10)  neglects  any  effects  of  eccentricity  in  the 
circularity  of  the  cylinder  and  of  lateral-torsional  stiffener  buckling  on  the 
compressive  stress  in  the  rings. 


FORMULATION  OF  OPTIMIZATION  PROBLEM 


As  indicated  in  the  Introduction,  the  optimization  problem  takes  the  follow¬ 
ing  three  alternative  forms: 

I.  Find  the  minimum-weight  design  of  the  stiffened  cylindrical  shell  subjected 
to  the  applied  loads  described  previously  and  constrained  against  overall 
(system)  buckling,  panel  (inter-ring)  buckling,  T  web  and/or  flange  buck¬ 
ling,  panel  and/or  ring  yield.  The  shell  is  also  required  to  possess  a 
lowest  natural  frequency  (in  vacuo)  greater  than  a  specified  minimum. 

II.  Find  the  structural  design  which  maximizes  the  separation  between  the  low¬ 
est  two  natural  frequencies  of  vibration  (in  vacuo)  for  the  stiffened  cylin¬ 
drical  shell  subjected  to  the  applied  loads  described  previously  and  con¬ 
strained  against  overall  (system)  buckling,  panel  (inter-ring)  buckling, 

T  web  and/or  flange  buckling,  panel  and/or  ring  yield.  The  shell  is  also 
required  to  possess  a  lowest  natural  frequency  greater  than  a  specified 
minimum,  and  a  structural  weight  less  than  a  prescribed  maximum. 

III.  Find  the  structural  design  which  maximizes  the  separation  between  the 
lowest  two  natural  frequencies  of  vibration  (in  vacuo)  which  have  primar¬ 
ily  axial  content  for  the  stiffened  cylindrical  shell  subjected  to  the 
applied  loads  described  previously  and  constrained  against  overall  (sys¬ 
tem)  buckling,  panel  (inter-ring)  buckling,  T  web  and/or  flange  buckling, 
panel  and/or  ring  yield.  The  shell  is  also  required  to  possess  a  lowest 
natural  frequency  greater  than  a  specified  minimum  and  a  structural  weight 
less  than  a  prescribed  maximum. 

Each  of  these  problems  has  the  form  of  an  inequality-constrained  optimiza¬ 
tion  problem,  which  may  be  solved  by  any  of  a  number  of  mathematical  programming 
algorithms.  The  particular  method  of  solution  chosen  for  this  work  is  the 
Sequential  Unconstrained  Minimization  Technique  (SUMT)  in  which  the  optimum 
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structural  design  problem  (in  either  form  I,  II,  or  III)  is  converted  into  a 

sequence  of  unconstrained  problems  [6]  by  means  of  so-called  "penalty  functions." 

In  this  technique,  it  is  required  to  find  a  vector  x  =  {x^x^. . . .x^}  = 

{tst(j)d(j)^xdftf ),  the  components  of  which  are  the  design  variables,  such  that  a 

specified  function  of  these  variables,  F(x),  is  extremized  while  satisfying  a 

set  of  constraints  of  form  g^(x)  £  0,  k  =  l,...,s. 

The  "objective  function,"  F(x),  takes  one  of  three  forms,  depending  on 

the  optimization  problem  being  considered.  For  Problem  I,  weight  minimization, 

27rp 

F (x)  =  W  =  - 2-s -  <RLt  +  (L/A  )[(R  -e  )t  d  +  (R  -e  )t  d  ]}  (12) 

TT (R+t  /2)2LP  S  X  *  *  41  f  £  f 

s  w 

where  parameters  e^  and  e^  are  defined  in  Appendix  I.  In  Eq.  (12)  F(x)  has 
been  normalized  by  dividing  shell  mass  by  the  mass  of  displaced  water. 

For  Problem  II,  the  case  of  separation  of  lowest  natural  frequencies,  an 
inverse  formulation  is  used.  The  separation  is  maximized  by  minimizing  the 
inverse  of  the  separation.  For  the  modes  being  considered,  an  ordered  list 


is  made  giving  U)  <  etc.  The  objective  function  is  then 


where  is  a  normalization  factor  taken  as  the  initial  frequency  separation. 
To  separate  frequencies  having  primarily  a  longitudinal  deformation  content 


(Problem  III)  it  is  necessary  to  examine  and  order  frequencies  in  each  primarily 
axial  mode,  i.e.,  the  ones  having  both  and  A^  smaller  than  A^.  The  objec¬ 
tive  function  is  then  given  by 


F(x) 


(14) 


where  the  subscript  L  has  been  added  to  denote  the  longitudinal  character  of 


the  shell  vibration. 
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The  number  of  design  variables,  q,  is  a  maximum  of  six  in  this  study,  but 
may  be  less  if  certain  of  the  design  variables  are  fixed.  Also  upper  and  lower 
] imits  IL  and  L^,  respectively,  are  assumed  specified  for  each  variable  x^. 

These  upper  and  lower  limit  constraints,  respectively,  are  written  in  the  normal¬ 
ized  form 


gi(x)  =  (U±  -  x±)/(Ui  -  Li)  £  0 

8q+i«  ■  <xi  -  Li>/(U1  '  Li>  S 


i  1,  . . . ,  q 


It  is  also  necessary  to  include  a  geometric  admissability  constraint  which 
serves  to  keep  the  frame  flanges  from  overlapping.  This  is  expressed  in  the 
normalized  form 

g2q+1(x)  -  1  -  df/*x  2  0  (17) 

The  behavioral  constraints  may  also  be  normalized.  The  panel  yield  con¬ 
straint  is  expressed  as 

S2q+2(S>  -  1  -  -  W  +  °x)1/2/0y  *  0  <18> 

The  frame  yield,  flange  buckling,  web  buckling  and  skin  buckling  constraints, 
respectively,  are  written  as 

g2q+3(2L>  -  1  -  l°cl/°y  *  0  (19) 

g2q+4W  ■  1  -  l°cl/0fcr  S  °  <20) 

82q+5(i>  =  1  -  |aci/0wcr  i  0  <2i; 

g2q+6«  -  1  -  */P„  *  0 

Finally,  to  prevent  gross  buckling  the  lowest  frequency,  oo^,  should  be 
greater  than  a  specified  minimum  frequency,  to  .  This  is  stated  in  the  form 

g2q+7(i>  =  (a,l  ’  a,min)/A2  *  0  (23> 


where  A^is  the  initial  value  of  to,  -  to  .  . 
L  1  min 
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For  the  minimum-weight  design  of  the  shell  these  are  all  the  constraints 

required,  but  in  order  to  separate  frequencies  two  other  constraints  must  be 

included.  It  is  conceivable  that  large  separations  could  be  achieved,  but 

possibly  only  at  the  expense  of  large  increases  in  weight.  It  is  therefore 

necessary  to  establish  an  upper  bound  to  the  mass  of  the  structure,  W 

max 

This  constraint  is  expressed  by 


g2q+8c*>  =  i 


W/W  >  0. 
max 


(24) 


The  final  constraint  is  not  required  in  the  definition  of  the  frequency 
separation  problems,  but  is  convenient  for  numerical  solution  when  penalty 
functions  are  utilized  in  the  Sequential  Unconstrained  Minimization  Technique. 
The  most  efficient  unconstrained  minimization  methods  require  the  computation 
of  gradients  of  the  objective  function  in  order  to  operate.  These  gradients 
should  be  smooth  and  continuous,  but  experience  shows  that,  as  03^  and  are 
separated,  and  03^  approach  a  common  value.  Eventually  and  03^  will 
switch,  with  the  result  that  the  mode  which  previously  represented  now 
represents  03^  and  vice  versa.  This  causes  discontinuity  in  the  gradient  of 
the  objective  function  and  subsequent  difficulties  with  the  numerical  algorithm. 
This  difficulty  can  be  overcome  simply  by  requiring  that  the  second  and  third 
frequencies  never  be  equal.  This  requirement  is  expressed  as  a  "singularity 
avoidance"  constraint  in  the  two  cases  of  frequency  separation  by 

82q+9(i)  =  “3  '  “2  >  0  (25) 

or 

g2q+9«  -  “l3  -  \2  >  0  • 

In  order  to  apply  the  methods  of  unconstrained  nonlinear  programming,  the 
basic  inequality-constrained  problem  is  recast  in  the  form  of  an  interior 
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penalty  function  [6],  the  solution  of  which  requires  finding  that  x  which 
minimizes 


$(x)  =  F(x)  +  r  .  (x) 
j=l  3 


(27a) 


where 


i/gjCx)  if  g.-(x)  -  £ 

l2e  -  gj(x)]/e2  if  gj(x)  <  G 


(27b) 


Pj (x)  is  the  so-called  extended  penalty  function  [7]  which  allows  the  use  of 
infeasible  designs  in  the  search  for  an  optimal  feasible  design.  The  quanti¬ 
ties  £  and  r  are  small  positive  scalars  and  s  is  the  number  of  constraints. 
Solving  this  unconstrained  minimization  problem  for  successively  smaller  values 
of  r  gives  a  series  of  designs  which  converge  to  a  local  or  global  optimum  of 
F(x).  This  procedure  is  referred  to  as  the  Sequential  Unconstrained  Minimi¬ 
zation  Technique  (SUMT).  Each  unconstrained  minimization  is  performed  by  the 
Davidon-Fletcher-Powell  method  [6].  In  order  to  apply  this  algorithm  it  is 
necessary  to  compute  the  gradients  of  the  objective  function  and  the  con¬ 
straints  with  respect  to  each  design  variable.  The  mathematical  complexity  of 
this  problem  makes  analytic  calculation  of  partial  derivatives  impractical 
and  it  is  found  convenient  to  use  a  forward  difference  technique  to  find  the 
gradients  numerically.  The  success  of  the  unconstrained  minimization  hinges 
on  the  ability  to  perform  an  accurate  unidirectional  minimization.  A  special 
hyperbolic  interpolation  formula  was  developed  [6]  for  use  with  the  SUMT 
method.  A  test  for  the  minimum  was  developed  and  incorporated  into  the  minimi¬ 
zation  algorithm  which  requires  that  a  measure  of  normality  between  the  direc¬ 
tion  vector  s_,  and  the  gradient,  V$,  be  less  than  0.001,  i.e., 

(s  •  V$)/(|s]  •  |V$|)  <  0.001  .  (28) 
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EXAMPLES 


Design  examples  given  in  [3]  were  re-evaluated  in  this  study  for  optimal 

performance  for  situations  I,  II  or  III  as  detailed  above.  In  all  of  these 

designs  the  preassigned  parameters  were  given  the  following  values: 

R  =  198  in.,  L  =  594  in. ,  E  =  30  x  106  lb/in.2,  Q  =  60  x  103  lb/in.2, 

-4  3 

p  =  7.75  x  10  slug/in.  ,  V  =  0.33  and  specified  pressure  associated  with 
depths  of  water  of  1000,  2000  and  3000  ft.  The  prescribed  minimum  natural 
(in  vacuo)  frequency  was  taken  as  12.0  Hz.  Except  for  the  dynamic  effects 
these  problems  are  quite  similar  to  those  in  [3],  The  structure  was  designed 
initially  for  minimum  weight  (problem  I)  for  the  three  different  operating 
depths.  The  results,  given  in  Table  1,  indicate  normalized  minimum  weights 
of  0.13317,  0.22295,  and  0.31922  for  operating  depths  of  1000,  2000  and 
3000  ft.,  respectively,  values  somewhat  lower  than  reported  in  [3],  This 
occurrence  is  due  to  the  fact  that  in  this  study  the  ring  spacing  was 
represented  by  a  continuous  variable  and  web  and  flange  thicknesses  were 
included  as  independent  design  variables,  while  in  [3]  the  frame  spacing  was 
a  discrete  parameter  and  the  web  and  flange  thicknesses  were  linked  and  re¬ 
quired  to  be  not  less  than  1/18  of  the  web  and  flange  depths.  In  the  designs 
presented  herein  the  frame  webs  are  very  thin  and  are  critically  stressed, 
i.e.,  on  the  verge  of  buckling.  Consideration  of  lateral-torsional  buckling 
and  the  effects  of  hull  eccentricity  may  alter  this  condition,  although  these 
effects  were  excluded  in  the  present  work. 

One  important  benefit  in  obtaining  the  minimum-weight  designs  in  Table 
1,  is  that  they  serve  as  initial,  feasible  designs  for  design  problems  II  and 
III,  provided  the  same  minimum  frequency  constraint  is  employed  and  the  maximum 
allowable  structural  mass  is  greater  than  that  of  the  minimum-weight  (mass) 
design  for  the  static  load  conditions. 
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Table  1.  Design  Problem  I  -  Weight  Minimization 
(R  =  198  in.,  L  =  594  in.,  E  =  30  x  106  psi,  0 =  60,000  psi) 


Depth 

1000  ft. 

2000  ft. 

3000  ft. 

t  ,  in. 
s* 

1.2108 

2.4856 

3.5156 

t.,  in* 

<f> 

0.3765 

0.4207 

0.4543 

d,  ,  in* 

<P 

19.589 

19.284 

19.915 

&  ,  in. 

33.602 

51.528 

36.195 

df,  in. 

17.664 

14.999 

16.991 

tf,  in. 

0.4705 

0.4984 

0.5453 

OJ^,  Hz. 

28.12 

26.05 

26.64 

^2 ,  Hz. 

49.39 

30.04 

35.88 

0)^ ,  Hz . 

52.31 

51.98 

55.11 

Weight3 

(Normalized) 

0.13317 

0.22295 

0.31922 

Maximum  (upper  bound)  weights  (normalized)  are  0.15,  0.25 
and  0.35,  respectively. 
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The  results  of  design  problem  II,  optimization  for  maximum  frequency 
separation  are  given  in  Table  2  for  the  cases  of  the  same  three  preassigned 
operating  depths.  The  maximum  allowable  normalized  weight  was  taken  as 
approximately  10%  greater  than  that  for  the  minimum-weight  design  for  the 
static  load  condition,  the  values  being  0.15,  0.25  and  0.35,  respectively, 
for  the  operating  depths  of  1000,  2000  and  3000  ft.  For  the  three  depth 
requirements  the  frequency  separation  was  increased  from  the  minimum-weight 
design  values  of  21.27,  3.993  and  9.2418  Hz  to  23.59,  24.28  and  25.33  Hz. 

Of  interest  is  the  fact  that  the  first  and  third  solutions  (for  1000  ft. 
and  3000  ft.  operating  depths)  gave  frequencies  =  ^3*  i.e.,  nearly  identi¬ 
cal  second  and  third  frequencies,  and  that  in  these  two  cases  the  maximum 
weight  constraint  was  less  than  critical.  It  is  thus  apparent  that  the 
major  portion  of  the  frequency  separation  is  obtained  by  bringing  the  second 
frequency  up  to  the  third  frequency,  and  that  little  or  nothing  is  to  be  gained 
by  adding  more  material  to  the  structure  after  this  is  accomplished.  In  the 
example  with  2000  ft.  depth  the  second  and  third  frequencies  were  unable  to 
completely  approach  each  other,  before  violating  the  maximum  weight  constraint, 
which  became  critical  in  this  case.  It  may  be  noted  that  for  these  designs 
the  frame  webs  are  very  thin  and  the  frame  flanges  are  relatively  thick.  It 
thus  seems  that  the  frequency  separation  has  been  achieved  by  making  the 
moment  of  inertia  of  the  frames  as  large  as  possible. 

Three  problems  of  category  III  were  designed  to  find  the  maximum  fre¬ 
quency  separation  in  primarily  longitudinal  modes  of  vibration  for  the  three 
different  operating  depths.  Results  are  given  in  Table  3.  For  the  cases 
which  have  operating  depths  of  1000  ft.  and  2000  ft.  the  algorithm  became 
entrapped  in  a  singularity  in  the  design  space.  The  problem  run  at  3000  ft. 
depth  was  unable  to  reach  the  singularity  because  the  maximum  weight 


17 


(R 


Table  2.  Design  Problem  II  -  Frequency  Separation 
=  198  in.,  L  =  594  in.,  E  =  30  x  10^  psi,  a =  60,000  psi) 


Depth 

1000  ft. 

2000  ft. 

3000  ft. 

t  ,  in. 

1.2216 

2.4717 

3.3986 

s’ 

t.,  in* 

<P 

0.3950 

0.3884 

0.5168 

dA,  in. 

20.722 

19.674 

23.764 

i  ,  in. 

33.853 

51.417 

36.065 

df,  in. 

17.551 

15.075 

16.864 

tf,  in. 

0.4653 

1.8638 

0.9485 

U)1,  Hz. 

28.3711 

29.4905 

29.7225 

u)2,  Hz. 

51.9638 

53.7674 

55.0508 

oo^,  Hz. 

51.9640 

54.1466 

55.0512 

Weight3 

0.1351 

0.25 

0.32928 

(Normalized) 

(u^-o^)  Hz. 

23.5928 

24.2769 

25.3283 

Maximum  (upper  bound)  weights  (normalized)  are  0.15,  0.25 
and  0.35,  respectively. 
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Table  3.  Design  Problem  III  -  Longitudinal  Frequency  Separation 

(R  =  198  in.,  L  =  594  in.,  E  =  30  x  10^  psi,  a ^  =  60,000  psi) 


Depth 

1000  ft. 

2000  ft. 

3000  ft. 

V  ln- 

1.5975 

2.6720 

3.9377 

y  in- 

0.2500 

0.4171 

0.42116 

y in- 

11.603 

17.384 

19.817 

i  ,  in. 

38.672 

52.232 

36.143 

df,  in. 

9.1454 

13.232 

16.907 

t^,  in. 

0.3826 

0.5248 

0.4927 

03^,  Hz. 

163.0006 

163.0640 

138.2134 

0J2,  Hz. 

192.2375 

192.7389 

162.4559 

Hz. 

222.6656 

222.7126 

221.6170 

Weight3 

0.14215 

0.23559 

0.34907 

(Normalized) 

(a>L  )  Hz. 

2  X 

29.2369 

29.6749 

24.2425 

Maximum  (upper  bound)  weights  (normalized)  are  0.15,  0.25, 
and  0.35,  respectively. 
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constraint  had  become  active.  A  singularity  avoidance  constraint  could  have 
been  developed  to  enable  the  algorithm  to  proceed,  however  this  was  not  done 
because  in  order  for  the  optimization  problem  to  have  a  significant  physical 
purpose,  a  better  definition  of  what  actually  constitutes  a  longitudinal  mode 
is  needed.  As  may  be  seen  in  Table  3  for  the  3000  ft.  case,  the  mode  having 
one  longitudinal  wave  and  no  circumference  waves  (m  =  1,  n  =  0)  has  a  "longi¬ 
tudinal"  frequency  of  138.21  Hz,  because  the  v-component  of  the  associated 
eigenvector  is  zero  and  the  w-component  is  less  than  1.0  (0.977).  As  the 
algorithm  separates  the  frequencies,  the  eigenvector  associated  with  this 
frequency  changes  form  with  the  w-component  increasing  to  a  value  slightly 
greater  than  1.0.  As  this  occurs,  a  second  eigenvector  of  this  mode’s 
vibratory  class  (m  =  1,  n  =  0)  becomes  longitudinal  in  form,  according  to 
the  definition  currently  in  use.  The  result  is  the  situation  encountered 
for  the  1000  ft.  and  2000  ft.  cases,  where  the  m  =  1,  n  =  0  mode  has  the 
second  lowest  frequency.  The  algorithm  becomes  entrapped  at  a  frequency 
separation  of  29  Hz  as  two  eigenvectors  for  the  m  =  1,  n  =  0  mode  switch 
back  and  forth,  having  w-components  both  approximately  equal  to  1.0.  Since 
the  m  =  1,  n  =  0  mode  has  two  frequencies  of  vibration  with  almost  the  same 
mode  shape  it  is  not  realistic  to  call  either  one  the  unique  longitudinal 
frequency  for  that  mode. 

This  effect  was  only  recently  encountered  and  must  be  given  additional 
study.  When  successfully  resolved,  the  result  will  be  a  capability  for 
generating  optimum  solutions  to  design  problems  I,  II  and  III.  For  the 
future,  research  should  be  directed  toward  incorporating  more  sophisticated 
structural  models,  with  more  design  flexibility,  in  anticipation  of  obtaining 
even  more  efficient  structural  configurations  from  both  the  standpoint  of 
minimum-weight  and  maximum  frequency  separation. 
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APPENDIX  I:  DYNAMIC  ANALYSIS 

The  forces  on  the  cylindrical  shell  are  as  shown  in  Fig.  2.  The  forces 
are  expressed  in  terms  of  displacements  as: 

Nx  =  Hu’  +  (Hv/R) v*  -  (Hv/R)w  +  (D/R)w" 

N^  =  HVu*  +  [ (H+HC+HF) /R] v*  -  [H+HC (1+e^/R)  +  HF(l+ef /R) ]w/R 

-  [D/R+HC (e^+p^/R)  +  HF(ef+p^/R)]^2  w** 

R 

=  (S/R)u*  +  Sv’  +  (K/R2)wf* 

N.  =  (S/R)u*  +  Sv’  -  (K/R2)w’ * 

(px 

Mx  =  -  (D/R)u'  -  (Dv/R2) v*  -  Dw"  -  (Dv/R2)w** 

=  (HCe^/R  +  HFef/R)v*  -  [D/R  +  HC(e^+p2/R)  +  HF(ef+p2/R) ]w/R 

-  Dvw"  -  [D  +  HC(p2+a3/R)  +  HF(p2+a2/R)w**/R 
Mx<()  =  -  (2K/R)  (v'+w' *) 

=  (K/R2) u*  -  (K/R)v'  -  (2K+Q)w’*/R 


(  )*  = 


3(  ) 


(),  =in 


8<j>  ’  v  '  "  3x 


The  matrices  [K] ,  [K  ]  and  [M]  of  equation  (6)  are: 

(j 

k 

[K]  =  k 
k 


11 

k12 

k13 

12 

k22 

k23 

13 

k23 

k33 

2  2  2 
kll  =  -  ^  “  (S/R  )n 


k12  =  “  (Hv+S) An/R 

k13  =  "  tHvA+DA3  -  (K/R2) An2] /R 

k22  =  "  (s+2k/r2>a2  -  [H+HC (1-e^/R)  +  HF(l-ef /R) ]n2/R2 
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k23  =  -  (3K+Dv)A2n/R2 

+  [D/R2-H+HC (p2/R2-l)  +  HF(p2/R2~l) ]n/R2 

-  [HC  (a2/R2-e^)  +  HF(aj|/R2-ef)]n3/R3 
k33  =  -  [H+HC (1+e^/R)  +  HF(l+ef/R) ]/R2 

+  2[D/R+HC (e^+p2/R)  +  HF(ef+p2/R)]n2/R3 

-  DA4  -  [ 2DV+4K+Q] (An/R) 2 

-  [D4-HC  (p2+a3/R)  +  HF(p2-Kx3/R)](n/R)4 


kGll 

0 

kG13 

=  0 

kG22 

kG23 

_kGl3 

kG23 

kG33 

CM 

+ 

CM 

c 

II 

/  2)  /R 

kG23  = 

n/R 

=  -X 

lc  = 

G33 

-  (n2+A2/2)/R 

kG22  =  "  (n2+x2/2)/R 

mll  0  m13 

M  =  0  m22  m23 

m13  m23  m33 

m  m> 

rall  =  ts  +  A(1-e4>/R)  +  B(l-ef/R) 
ml3  =  ttg/(12R)  +  A(-e^+p2/R)  +  B(-ef+p2/R) ]A 
t2 

m22  =  ts(1  +  ^2  }  +  A(l-3e(})/R+3p2/R2-a3/R3) 

+  B(l-3ef/R+3p2/R2-a3/R3) 
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m23  =  [tf/(6R2)  +  A(-e^/R+2p2/R2-a3/R3) 

+  B(-ef/R+2p2/R2-a3/R3)]n 

m33  =  ts  +  Cx2+n2/R2) tg/12  +  Atl-e^/R  +  (A2+n2/R2) (p2-a3/R) ] 
+  B[l-ef/R+(A2+n2/R2)(p2-a3/R)] 


The  section  properties  are  defined  as  follows: 


H  =  Et  /(1-v  ) 
s 


D  =  Et3/[12(1-V2)] 
s 


A  =  t.d  Jl 

cj)  <f>' X 

HC  =  EA 
G  =  E/ [2 (1+V) ] 


B  =  t.d./fl, 
f  f  x 


HF  =  EB 
t  =  t  +  2d 


S  =  Gt 


K  =  Gt  /12 
s 


Q  =  G(J^+ Jf)Aa 

J,  =  c,d ,t3 
9  9  9  9 


Jf  =  CfdfC£ 


c ,  =  -  0.285  e 
9 


-(0.49  A 


Cf  =  -  0.285  e 


-(0.49 


V  +  0.316 
tf/  +  0.316 


e ,  =  y(d,+t  ) 
9  2  cp  s 


ef  =  2(df+tt) 


2  L2  ,  L  ,12 

9  3  (p  2  cp  s  4s 

2  1.3 , 1  ,12 

f  3  f  2ft  4t 


3  1j3  1  2  3  2.  ^  13 

Ot .  =  yd  ,  +  — t  d  .  4-  -yt  d  ,  +  yt 

9  4  cp  2  s  <p  8  s  9  8s 

3  1J3  ,  1  .232,  ,13 

af  4^  +  2ttdf  +  8ttdf  +  8ts 
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APPENDIX  II:  STATIC  STRENGTH  ANALYSIS 
The  critical  compressive  stresses  in  the  skin  are  assumed  to  occur  on  the 
outer  surface  at  mid  panel.  The  stresses  there  are 
°4>  -  -  (ER/ts)[l+r(Hn+VHE)] 

°X  -  -  <PR/V  [T  +  ne] 


where  the  various  parameters  are  given  as: 
-PR/tg  =  hoop  stress 

R  =  radius  to  mid  plane  of  shell 

T  =  [ 1— v/ 2— B] / (1+3) 


The  hoop  compressive  stress  in  the  frame  is 
ac  =  PVR/ A 

where  V  is  a  load  factor  given  as: 

V  =  t^[l+(l-v/2) 3/B] / (1+3) 
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